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ABSTRACT 

A typical galaxy is thought to contain tens of millions of stellar-mass black 
holes, the collapsed remnants of once massive stars, and a single nuclear su- 
permassive black hole. Both classes of black holes accrete gas from their envi- 
ronments. The accreting gas forms a flattened orbiting structure known as an 
accretion disk. During the past several years, it has become possible to obtain 
measurements of the spins of the two classes of black holes by modeling the X-ray 
emission from their accretion disks. Two methods are employed, both of which 
depend upon identifying the inner radius of the accretion disk with the innermost 
stable circular orbit (ISCO), whose radius depends only on the mass and spin of 
the black hole. In the Fe Ka method, which applies to both classes of black holes, 
one models the profile of the relativistically-broadened iron line with a special 
focus on the gravitationally redshifted red wing of the line. In the continuum- 
fitting method, which has so far only been applied to stellar-mass black holes, 
one models the thermal X-ray continuum spectrum of the accretion disk. We 
discuss both methods, with a strong emphasis on the continuum-fitting method 
and its application to stellar-mass black holes. Spin results for eight stellar-mass 
black holes are summarized. These data are used to argue that the high spins of 
at least some of these black holes are natal, and that the presence or absence of 
relativistic jets in accreting black holes is not entirely determined by the spin of 
the black hole. 

1. Introduction 

Our focus throughout is on the two principal methods for measuring the spins of accret- 
ing black holes: modeling the thermal continuum X-ray spectrum and modeling the profile of 
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the relativistically-broadened Fe Ka line. The continuum-fitting (CF) method, has thus far 
only been applied to stellar-mass black holes in X-ray binaries, whereas the Fe Ka method 
has been applied to both stellar- mass (M ~ 10 Mq) and supermassive (M ~ 10® — 10^" M^ 



black hole^j. This paper is chiefiy focused on measuring the spins of stellar-mass black holes 
via the CF method because we have been deeply engaged in this work during the past sev- 
eral years. In Section 6, we secondarily discuss the Fe Ka method, which is very important 
because it is the primary approach to measuring the spins of supermassive black holes. 

We note that spin estimates have also been obtained by modeling the high-frequency 



X-ray oscillations (100-4 50 Hz) observed for several stellar-mass black holes (IWagoner et al. 



2001 



Torok et al.l 120051 ). At present, this method is not providing dependable results be- 



cause the correct model of these o scillations is not known. X-ray polarimetry is anothe r 
potential avenue for measuring spin ( iDovciak et al.ll2008l : iLi et al.ll2009l : ISchmoU et al.ll2009l ). 
which may be r ealized soon with t he 2014 launch of the Gravity and Extreme Magnetism 
Small Explorer ( ISwank et al.ll2008l ). Meanwhile, several other method s of measuring spin 



have been proposed or applied (e.g.. lTakahashill2004l : iBarai et al.l |2004J : iHuang et al. 
Suleimanov et al. 2008 : Shcherbakov fc Huang! 2011 ). 



2007 



2. Stellar-Mass Black Holes in X-ray Binaries 

Observations in 1972 of the X-ray binary Cygnus X-1 provided the first strong evidence 
that black holes exist. Today, a total of 23 such X-ray binary systems are known that contain 
a compact object too massiv e to be a neutron star or a degenerate star of any kind (i.e.. 



M > 3 Mq] lOzel et al.ll2010l ). These compact objects, which have typical masses of lOM©, 
are referred to as black holes. Their host systems are mass-exchange binaries containing 
a nondegenerate star that supplies gas to the black hole via a stellar wind or via Roche- 
lobe overfiow in a stream that emanates from the inner Lagrangian point. The mass-donor 
star in the Roche-lobe overfiow systems is typically a low mass (M ~ 1 Mq) sun-like star, 
and the X-ray source is transient, alternating between yearlong outbursts (Lip^x ~ -^Edd = 
1.3x lO^^M/lOM© erg s~^) and years or decades of quiescence {L ~ 10^^LEdd)D- The wind- 
fed X-ray sources, on the other hand, are fueled by massive hot stars (M > IOMq) and 
are persistently luminous. A schematic sketch to scale of 21 of these systems is shown in 
Figure 1: The four at top are persistent systems and t he 17 at bottom are the transie nts. 



For a review of the properties of black hole binaries, see iRemillard fc McClintockl (120061 ). 



H Mq = 1 solar mass = 2.0 x lO^^ g. 

^^Edd is the critical Eddington luminosity above which radiation pressure exceeds gravity. 
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3. The Continuum-Fitting Method 



A definite prediction of relativity theory is the existence of an innermost stable circular 
orbit (ISCO) for a particle orbiting a black hole. This inherently relativistic effect has a 
major impact on the structure of an accretion disk. At radii R > -Risco (the radius of 
the ISCO), accreting gas moves on nearly circular orbits and slowly spirals in toward the 
black hole. At the ISCO, however, the dynamics change suddenly and the gas, finding no 
more stable circular orbits, plunges into the hole. In the continuum-fitting (CF) method, 
one identifies the inner edge of the accretion disk with the ISCO (see Sees. 4 and 5 for 
supporting evidence) and estimates -Risco by fitting the X-ray continuum spectrum. Since 
the dimensionless radius risco = Risco /(GM/c^) is solely a monotonic function of the black 
hole spin parameter a* (Fig. 2)|f|, knowing its value allows one immediately to infer the value 
of a*. We note that the truncation of the disk at the ISCO is also a crucial assumption of 
the Fe Ka method of measuring spin ( iReynolds &: Fabianll2008[ ). 



Before describing the CF methodology, we stress that for this technique to succeed it 
is essential to have accurate measurements of the distance D to the source, the inclination 
i of the accretion disk, and the mass M of the black hole (for reasons discussed below). 
The methodologies for measuring D, i and M are firmly established. Therefore, rather than 
digressing to discuss how these me asurements are rnade, we refer the interested read er to 
some recent papers on the subject (lOrosz et al.l 120071 . l2009l . l2011t ICantrell et al.l 120101 ). 



The gaseous matter fiowing from the companion star to the black hole has appreciable 
angular momentum as a consequence of the binary orbital motion. As the gas fiows, viscous 
forces cause it to spread out into an orbiting structure known as an accretion disk. The gas 
fiowing into the outer disk spirals slowly inward on Keplerian orbits on a time scale of weeks, 
reaching a typical temperature near the ISCO of kT ~ 1 keV. Because the accretion disk is 
of fundamental importance to the measurement of black hole spin, we now describe in some 
detail the thin-disk model we employ. 

The model we use is that described in iNovikov fc Thornel Jl973[ hereafter NT) , which i s 
a relativisitic generalization of a Newtonian model developed by lShakura fc SunyaevI ( 119731 ). 
The NT model describes an axisymmetric radiatively-efficient accretion fiow which, for a 
given black hole mass M, mass accretion rate M and black hole spin parameter a*, gives 
a precise prediction for the local radiative fiux F{R) emitted at each radius R of the disk. 
Moreover, the accreting gas is optically thick and the emission is thermal and blackbody-like, 



■^We express black hole spin in the customary way as the dimensionless quantity a* = cJ/GM^ with 
I a* I < 1, where M and J are respectively the black hole mass and angular momentum. 
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making it straightforward to compute the spectrum of the emission. Most importantly, the 
inner edge of the disk is located at the ISCO. Therefore, from the measurement of -Risco; 
and if we know the mass M of the black hole, we can immediately obtain a* (Fig. 2). This is 
th e principle behind t he CF method of estimating black hole spin, which was first described 



bv lZhang et all (Il997f ). 



Before discussing how to measure -Risco of a disk, we remind the reader how one mea- 
sures the radius i?* of a star. Given the distance D to the star, the radiation flux Fobs 
received from the star, and the temperature T of its continuum radiation, the luminosity of 
the star is given by 

= AttD^Foi^s = 47^i^>T^ (1) 

where a is the Stefan-Boltzmann constant. From Fobs and T, we immediately obtain the 
quantity n{R^/Dy, which is the solid angle subtended by the star. From this and the 
distance D, we immediately obtain R^. For accurate results we must allow for limb darkening 
and other non-blackbody effects in the stellar emission by computing a stellar atmosphere 
model, but this is a detail. 

The same principle applies to an accretion disk, but with some differences. First, since 
the flux F{R) emitted locally by the disk varies with radius R, the radiation temperature 
T{R) also varies with R. But the precise variation is known (if we assume the NT disk 
model), so it is easily incorporated into the model. Second, since the bulk of the emission is 
from the very inner regions of the disk, the effective area of the radiating surface is directly 
proportional to the square of the disk inner radius, Agg = CRf^^Q, where the constant C is 
known. Third, the observed flux Fobs depends not only on the luminosity and the distance, 
but also on the inclination i of the disk to the line of sight. Allowing for these differences, 
one can write a relation for the disk problem similar in spirit to eq. (1), i.e., given Fobs and 
a characteristic T (from X-ray observations), one obtains the solid angle subtended by the 
ISCO: TTcosi (-Risco/-D)^. If we further know i and D, we obtain -Risco! and if we also know 
M, we obtain a* (Fig. 2). This is the basic idea of the CF method. 

There are three main issues that must be dealt with before applying the method: (1) 
One must carefully trace rays from the disk to the observer in the Kerr metric of the rotating 
black hole in order to compute accurately the observed flux and spec trum. To this end, we 
have developed an accretion disk model called kerrb^ ( Li et al. for fitting X-ray dat 



a. 



KERRBB assumes the NT model of the disk and carries out all the necessary ray-tracing to 



''This model name and those that follow designate publicly-available 
comprise a suite of X-ray data analysis software known as XSPEC 



TOgrams that 
Arnaudl Il996t 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.htmll. 
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relate disk properties to observables. (2) One must have an accurate model for computing the 
spectral hardening factor / = T/Tcg, where T is the temperature of the radiation at a given 
radius and T^s is the effective temperature at the same radius defined by F{R) = aT^g{R). 
This correction for non-blackbody effects is important at the high temperatures typically 
found in b lack hole disks. To carr y out this correcti o n we use the advanced disk atmosphere 
models of iDavis et al.l (120051 ) and lDavis fc Hubenyl ( l2006l ). (3) Most importantly, the inner 
accretion disk must be well described by the standard geometrically-thin, optically-thick NT 
disk model that we employ. To ensu re this, we restrict ou r attention strictly to observations 



with a strong thermal comp onent (ISteiner et a. 



30% of the Eddington limit flMcChntock et al. 



2009aJ) and with disk luminosities below 



20061 ). In the two sections that follow, we 



present observational and theoretical evidence that at these luminosities the NT model is 
quite accurate. 



For a full description of the mechanics of the CF method we refer the reader to Section 4 



m 



McClintock et al. 



( I2OO6I ). In brief, we fit the broadband X-ray continuum spectrum in 
conjunction with other components as needed, principally a Compton tail component that 
i n more recent work is described using an empirical model of Comptonization called SIMPL 
( ISteiner et al.l l2009bl ). Th e accretion-dis k component, which is key for the CF method, 
is modeled using kerrbb ( iLi et al.ll2005[ ). which includes all relativistic effects within the 
context of the NT model, and also incorporates the advanced treatment of spectral hardening 
mentioned above. It furthermore includes self-irradiation of the disk ( "returning radiation" ) 
and the effects of limb darkening. The key fit parameters returned are the black hole spin 
parameter a^, and the mass accretion rate M. 



4. Truncation of the Accretion Disk at the ISCO: Observational Evidence 

A crucial assumption that underlies both the CF and Fe Ka approaches to measuring 
spin is that the accretion disk is quite sharply truncated at the ISCO. This assumption, 
which is fundamental to the NT model, is clearly valid if one considers only geodesic forces 
in the midplane. However, there are strong magnetohydrodynamic (MHD) forces acting 
in black hole accretion disks, and it is therefore unclear a priori that the disk terminates 
sharply at the ISCO. In this section, we present the observational evidence that there exists 
a fixed inner-disk radius in black hole binaries and, in the following section, we discuss the 
theoretical evidence for identifying this radius with rigco- 



There is a long history of evidence suggesting that fitting the X-ray continuum is a 
promising approach to measuring black hole spin. This history begi ns in the mid-1980s 
with the application of a simple nonrelativistic multicolor disk model ( iMitsuda et al.lll984j : 
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Makishima et al.lll986[ ). now known as diskbb, which returns the color te mperature at 
the inner-disk radius Ri^. In their review paper on black hole binaries, iTanaka fc Lewin 



( 119951 ) summarize examples of the steady decay (by factors of 10-100) of the thermal flux of 
transient sources during which Ri^ remains quite constant (see their Fig. 3.14). They remark 
that t he constancy of R, n suggests that this fit parameter is related to the radius of the 
ISCO. IZhang et al.l ( 119971 ) then outlined how, using a relativistic disk model and corrections 
for the effects of radiative transfer, the fixed inner disk radius provides an observational 
basis to infer black hole spin. More recently, the evidence for a constant inner radius in 
the thermal state has been presented for a number of sources via plots showing that the 
bolometric luminosity of the thermal component is approximately proport i onal to T;t (e.g. 



Kubota et al. 



2001 



McClintock et al. 



Ku 



2009 



Dunn et al.ll2010[ ) 



3ota fc Makishimal 120041: iGierhhski fc Pond 120041 : lAbe et al.l 12005 



A recent study of the persistent sour ce LMC X-3 presents the most compelling evidence 
to date for a constant inner-disk radius ( ISteiner et al.ll2010al ). We analyzed many spectra 
collected during eight X-ray missions that span 26 years. As shown in Figure 3, for a selected 
sample of hundreds of spectra obtained using the Rossi X-ray Timing Explorer (RXTE), we 
find that to within ^ 2 percent the inner radius of the accretion disk is constant over 
time and unaffected by the gross variability of the source (top panel). Meanwhile, even 
considering an ensemble of eight X-ray missions, we find consistent values of the radius to 
within ^ 4 percent. These results provide compelling evidence for the existence of a fixed 
inner-disk radius and establish a firm foundation for the measurement of black hole spin. 
The only reasonable inference is that this radius is closely associated with the radius of the 
ISCO, as we show to be the case in the following section. 



5. Truncation of the Accretion Disk at the ISCO: Theoretical Evidence 

The relativistic NT model on which the CF method is currently built assumes, as does its 
predecessor the Newtonian model of Shakura & Sunyaev (1973), that the accretion disk under 
consideration is geometrically thin. That is, the model assumes that the vertical thickness 
H at any radius R satisfies H <^ R. Assuming in addition that the disk is axisymmetric 
and in steady state, the model derives a number of relations which follow directly from basic 
conservation laws. One of the powerful results of this analysis is a formula for the disk flux 
profile F{R) which depends only on the mass M and spin a* of the black hole and the mass 
accretion rate M, but is independent of messy details such as the viscosity of the accreting 
gas. It is the existence of this robust result for F{R) that enables the CF method to work 
so well. 
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There is, however, one unproven assumption in the NT model which is incorporated via 
a boundary condition: The model assumes that the shear stress (which drives the accretion at 
larger radii) vanishes at the ISCO. This "zero-torque" assumption is intuitively reasonable 
(since the gas switches to a rapidly plunging state once it crosses the ISCO, why should 
there be a st ress at the t ransit ion r adius?), but as NT theniselve s realized, it is ultimately an 
assumption. iPaczyhskil (120001 ) and lAfshordi fc Paczyhskil (120031 ) argued that deviations from 
the NT model decrease monotonically with decreasing disk thickness and that thin disks with 
H/R <ti 1 should be very well described by the model. Their argument, which was based 
on a h ydrodynamical descri ption of the disk, was confirmed by detailed calculations carried 
out by IShafee et al.l ( l2008bl ). However, this leaves open the question of whether magnetized 
disks might deviate substantially from the NT model. In their paper, NT explicitly mention 
that magnetized disks could very well violate the zero-torque boundary condition. 

Arguments have been advanced to s uggest that a magnetized ac creting gas will indeed 



have a non-zero shear stress at the ISCO (jKroliklll999l : lGammielll999l ). and that furthermore 



this stress could be so large that it may completely invalidate the NT model even in very 
thin disks. This is clearly an important question that strikes at the heart of the CF method. 



tic magnetohydrodynamic (GRMHD) simulations have explored this question ( 


Shafee et al. 


2008a: 


Revnolds & Fabian 


2008; 


Noble & Krohk 


2009; 


Noble et al. 


2010; 


Penna et al. 


2010 





The conclusion of these authors is that the shear stress and the luminosity of the simulated 
disks do differ from the NT model, but perhaps not by a large amount. 



Figure 4, taken from iKulkarni et al.l ( 120101 ). shows the disk luminosity distribution 



Penna et al. 


c 


2010) 


Kulkarni et al 


. for 



0.045-0.08 (see Penna et al. and 



clear from the figure, the numerical models follow the NT model reasonably well, although 
they do deviate from it. Two kinds of deviation are seen. First, the numerical models pro- 
duce some radiation inside the ISCO, whereas the NT model predicts no radiation there. 
Second, the peak of the emission in the simulated disks is shifted inward relative to the peak 
in the NT model. Both of these effects cause the disk to appear to have a smaller ISCO 
radius. This in turn means that, if one were to infer the black hole spin by fitting this lu- 
minosity distribution (or the corresponding spectrum) using the NT model, one would infer 
an erroneously large value for the spin. This systematic error arises because the NT model 
is not a perfect description of the simulated disk. To the extent that the simulated disk is 
a closer match to a real accretion disk than is the NT model, this allows us to estimate the 
corresponding systematic error in our measurements of spin. 
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How serious is this systematic error? We answer this in three parts (see iKulkarni et aL 



20101 . for details). 



1. Each of the models shown in Figure 4 causes an error in the spin estimate that is 
smallest for a low disk inclination (face-on disk; see Fig. 1) and largest for a high disk incli- 
nation {i = 75°). For the latter (most unfavorable) case, the four models, which correspond 
to true spin values of a* = 0, 0.7, 0.9 and 0.98, give via the CF method spin values of 0.17, 
0.83, 0.936 and 0.991, respectively. 

2. A similar exercise can be carried out for disks with other thicknesses. It is found 
that the error in the spin estimate is larger for thicker disks and smaller for thinne r disks . 



Very roughly, the erro r appe ars t o scale as HI R. Thus , the conclusion of iPaczyhskil ( 120001 ) 



Afshordi fc Paczyhskil (j2003[ ) and IShafee et al.l ( l2008bl ). that deviations from the NT model 



vanish in the limit of vanishingly small disk thickness, appears to be valid also for MHD 
disks. 

3. The particular models shown in Figure 4 correspond to disk luminosities in the range 
L/L^dd ~ 0.4 — 0.8, based on their H/R. (It is difficult to be very quantitative since the 
mapping between luminosity and H/R is not known precisely.) Since the CF method is 
applied only to observations at L/L-^dd < 0.3, the systematic error due to inaccuracies in 
the theoretical model could be up to a factor of 2 smaller than the errors quoted in point 1 
above. 

Although the above results are based on numerical simulations that do not necessarily 
mimic real disks perfectly, we believe they still provide an estimate of the likely magnitude 
of the systematic error. The key point is that the level of systematic error we find is not 
serious at the current time. The observational errors considered in the following section are 
significantly larger in all cases. Stating this differently, while magnetized disks do behave as 
if their inner edges are shifted inward relative to the position of the ISCO (Fig. 4), the effect 
is quantitatively not serious for the disk luminosities (or disk thicknesses) at which the CF 
method is applied. 



6. Results of Continuum Fitting 

The spin results obtained to date for eight stellar-mass black holes are summarized 
in Table 1. The spins we find cover the full allowable range of prograde spins (Fig. 2) 
from a^, ^ (Schwarzschild) to a* ^ 1 (extreme Kerr). Interestingly, the spin values that 
have been obtained (Table 1), while spanning the full range of prograde spins, are all in 
the canonical physical range, namely |a*| < 1. This is an important result and not at all a 
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foregone conclusion: Given the hard external constraints on the dynamical model parameters 
(Z), i and M; Section 3), it is entirely possible that a black hole will be found that implies 
a spin beyond the reach of our current model, i.e., \a^\ > 1. Such a result could simply be 
caused by large systematic errors in D, i and M. Or, more interestingly, it could falsify our 
model, or even possibly point to new physics. This is why we are excited about improving 
our measurements of D, i and M for the near-extreme Kerr hole GRS 1915+105 (see Sec. 
6.2). 

Error estimates are primitive for the first four spin results published in 2006 (Table 1). 
In recent work on the other four black holes, the principal sources of observational error, 
as well as the uncertainties in the key model parameters (e.g., the viscosity parameter), 
have been treated in detail. In particular, in our most recent paper on XTE J1550-564, 
we exhaustively exp lored many different sources of error (see Table 3 and Appendix A in 
Steiner et al.ll2010bl ). The upshot of the work to date is that in every case the uncertainty 
in a^, is completely dominated by the errors in the three key dynamical parameters that 
we input when fitting the X-ray spectral data. As discussed in Section 3, these parameters 
are the distance D, the black hole mass M, and the inclination of the inner disk i (which 



we a ssume is aligned with the orbital angular momentum vector of the binary; iLi et al. 



20091 ). In order to determine the error in due to the combined uncertainties in D, M and 
z, we perform Monte Carlo simula tions assuming t hat these parameters are normally and 
independently distributed (e.g., see lGou et al.ll2009l ). Note that the errors introduced by our 
use of the NT model (Sec. 5), which are not consi dered here, are signifi cantly smaller than 
the observational errors (see Table 7 and Sec. 4 in lKulkarni et al.ll2010l ). 



We now discuss the results for four black holes in some detail. We first consider the 
persistent source M33 X-7. We then turn to GRS 1915+105, the prototype microquasar 
(IMirabel fc Rodrfgueg|l994l ). which hosts a near-extreme Kerr hole. Finally, we consider the 
microquasars A0620-00 and XTE J1550-564, contrasting their behavior with that of GRS 
1915+105. 



6.1. M33 X-7: The First Eclipsing Black Hole 



A long observation of the galaxy M33 with the Chandra X-r ay Observatory led to the 
discovery of a black hole that is eclipsed by its companion star (jPietsch et al.ll2006l ). We 
made a detailed follow-up dynamical study of the optical counterpart of M33 X-7, the first 
such study of a black hole binary beyond the environ s of the Milky Wa y. We determined a 
precise mass for the black hole, M = 15.65 + 1.45 Mq (jOrosz et al.ll2007l ). as well as the mass 
of its exceptional companion star (^ 70 Mq). As we discuss in Orosz et al., it is difficult 
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to understand the origin of this system - a massive black hole in a 3.5-day orbit, separated 
by only 42 solar radii from its supergiant companion (see Fig. 1). Recently, a consistent 
evolutionary model ha s been proposed that accounts for all the key properties of the system 
( iValsecchi et al.ll2010l ). It assumes that M33 X-7 started as a primary of ~ 95 Mq and a 
secondary of ~ 30 Mq, with an orbital period that is close to its present 3.5-day value. 

Using as input our precise values for the black hole mass and orbital inclination angle, 
and the well-established distance of M33, we fitted 15 Chandra and XMM- Newton X-ray 
spectra and obt a ined a precise value for the spin of the black hole primary, = 0.84 ± 0.05 
( iLiu et al.ll2008l . |2010| ). Remarkably, given that an (uncharged) astrophysical black hole is 



described by just its mass and spin, this result yields a comjglete description of an asteroid- 
size object at a distance of 2.74 million light-years (840 kp 



nple 



What is the origin of the spin of M33 X-7? Was the black hole born with its present 
spin, or was it torqued up gradually via the accretion flow supplied by its companion? In 
order to achieve a spin of a* = 0.84 via disk accretion, an initially non-spinning black hole 
must accrete 5.7M0 from its donor star (IKing fc Kolblll999l ) in becoming the M = 15.65M0 
black hole that we observe today. However, to transfer this much mass even in the case of 
Eddington-limited accretion requires > 17 million year^, whereas the massive companion 



star, and hence its host system, can not be older than about 2-3 million years (jOrosz et al. 



20071 ). Thus, it appears that the spin of M33 X-7 must be chiefly natal - i.e., the event horizon 
trapped much of the angular momentum of the collap sing stellar core - a cq i iclusion that 



has b een reached for two other stellar- mass black holes (iMcClintock et al.ll2006l : IShafee et al. 
20061 ) . (However, see lMoreno Mendez et al.ll2008l on hypercritical accretion). 



6.2. GRS 1915+105: A Near-Extreme Kerr Black Hole 



GRS 1915+105 has unique and striking properties that sharply disting uish it from the 50 



or so known binaries t hat are believed to contain a stellar-mass black hole (IMcClintock et al. 



20061 : lOzel et al.ll2010l ) . It is the most relia ble source of relativistic rad io jets in the Milky Way 
and is the prototype of the microquasars (IMirabel &: Rodrfguezlll994l ). It frequently displays 

extraordinary X-ray variability that is not mimicked by any other black hole system. The 

prop erties of its high-frequency X-ray oscillations are equally extraordinary (IRemillard fc McClintock 
20061 ). Among the 17 transient black hole systems, GRS 1915+105 is unique in having re- 



^1 parsec (pc) = 3.26 light-years. 

^-^^Edd = L^dd/vc^j where L^dd = 1.3xlO^^M/10Mf;i erg and the efficiency ij increases from 5.7% to 



13.3% as a» increases from to 0.84 ( Shapiro fc Teukolskvll 19831) . 
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mained active for more than a decade since its discovery in 1992. The system has an orbital 
period of 30.8 days, and it is the largest of the black hole binary systems (Fig. 1). 

The pc-scale radio jets of GRS 1915+105, with apparent velocities greater than the 
speed of light (superluminal motion), are the analogue of the kpc-scale jets that have long 
been observed for quasars. The source episodically ejects material at relativistic speeds, 
which can easily be tracked for w eeks at centimeter waveleng t hs as clouds of plas ma moving 
outward on the plane of the sky (IMirabel fc Rodrfgueg|l994l : iFender et al.lll999l ). Based on 
a kinematic model, the jet velocity for a plausible distance of D ~ 11 kpc is f j/c > 0.9, and 
the inclination of the jet from our line of sight is ij ~ 65°. 

Based on the analysis of X-ray spectral data for reasonable estimates of D, M and 
i, we discovered t hat GRS 1915+105 contains a near-extreme Kerr hole with > 0.98 



(iMcClintock et al.l 120061 ). However, the current estimates of both D and M are poor: The 
distance is uncertain by a factor of !^ 2 (Fig. 5b), and the mass is uncertain by ~ 30% 
(M = 14.4 ± 4.4 Mq; ICreiner et allboOll : iHarlaftis fc Greinei]l20o3 ). 



Remarkably, the extraordinarily high spin of GRS 1915+105 and other properties of the 
source allow one - for this source only - to place tight constraints on the allowable values 
of M and D, which must lie within the triangular region shown in Figure 5a. Our spin 
model constrains the black hole's mass and distance to lie to the right of the slanted line 
(99% confidence level) because the model implies values of a* > 1 to the left of the line. 
(Our model is only valid for a* < 0.999). Distances > 12 kpc are ruled out by the kinematic 
model of the radio jets (Fender et al. 1999). The lower bound on M is an estimate and 
is based on work in progress: We (Danny Steeghs et al.) are in the act of obtaining near- 
infrared spectroscopic data at ESO's VLT Observatory that we fully expect will improve the 
measurement of M by at least a factor of two. When this result is in hand, we will have M 
and D constrained to lie within the small shaded triangle, thereby constraining the distance 
to lie within the range 9.5-12 kpc. 

As noted above, the distance is also highly uncertain. Furthermore, all measurements 
to date are model-dependent estimates with large, systematic uncertainties that are difficult 
to assess (Fig. 56). With Mark Reid, we are in the act of obtaining a model- independent 
trigonometric distance with an uncertainty of 10% via a parallax measurement (the gold- 
standard method in astronomy) using the Very Long Baseline Array (VLBA), a worldwide 
array of radio telescopes. We have made successful observations at four epochs in 2008-2010 
and anticipate that observations at several additional epochs will be required to reach our 
goal. Two hypothetical and possible outcomes of these VLBA observations are indicated in 
Figure 5a: VLBA2 would confirm our spin model and VLBAl would rule it out. 
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6.3. A0620-00 and XTE J1550-564: Two Schwarzschild-Like Black Holes 



The host systems of these two black holes are quite small (Fig. 1). The optical companion 
in A0620-00, which has an orbital period of only 0.3 days, is a star somewhat cooler than 
the Sun with about half its size and mass. During its yearlong X-ray outburst in 1975-1976, 
this nearby transient {D ^ 1 kpc) became the brightest celestial X-ray source ever observed 
(apart from the Sun). For several days, the flux at Earth from this source was greater than 
the combined flux of all of the hundreds of other X-ray binaries in our galaxy. During this 
period, A0620-00 was also a bright transient radio source, w hich was observ e d with the early 
radio telescopes of the day. A reanalysis of these data by iKuulkers et al.l (119991 ) indicates 
that multiple jet ejections occurred. The authors find that, like GRS 1915+105, the radio 
source was extended on parsec scales, and they infer a relativistic expansion velocity. 

The cool companion star in XTE J1550-564 also has a mass about half that of the 
Sun, alth ough its radius is about twice as great (Fig. 1); the orbital period of the system is 
1.5 days (lOrosz et al.ll201ll ). During its principal 1998-1999 outburst cycle, this transient 
source produced one of the most remarkable X-ray flare events ever observed for a black 
hole binary. For ^ 1 day, the source intensity rose fourfold and the flux in the dominant 
nonthermal component of emission rose by the same factor. Th en, just as quickly, the source 
intensity dechned to its pre-outburst le vel (ISobczak et al.ll2000r). Four days later, small-scale 
superluminal radio jets were observed (IHannikainen et al.ll2009l ): their separation angle and 
relative velocity link them to the impulsive X-ray flare. The s ubsequent detection of pc-scale 
radio jets in 2000 led to the discovery of relativistic X-ray iets (ICorbel et al.ll2002l ). All of the 
available evidence strongly indicates that these pc-scale X-ray and radio jets are associated 
with the powerful X-ray flare. 

Using our recently-determined estimat es of D, M and i for A0620-00 and XTE J1550- 
564 (jCantrell et al.ll2010l : lOrosz et al.ll201ll ). we fltted the X-ray spectra of these black holes 
and determined their spins (Table 1). Figure 6 shows a pair of fltted spectra for the latter 
source. The spectrum in the left panel is completely dominated by the thermal component 
and is therefor e ideal for the deter r ninatiq n of spin. Now, however, using our improved 
methodologies (ISteiner et al.l l2009al . l2010bl ). we are able to obtain useful and consistent 
values of spin as well for spectra that have a strong Compton component of emission, like 
the one shown in the right panel of Figure 6. The origin of this component is widely attributed 
to Compton upscattering of the soft disk photons by coronal electrons (see Sec. 7). 

The CF spins of both A0620-00 and XTE J1550-564 are quite low: a, ^ 0.1 and 
a^, ^ 0.3, respectively. The corresponding nominal radii of their ISCOs are 5.7M and 5.0M, 
which differ only modestly from the Schwarzschild value of 6M. The low spins of these 
two microquasars challenge the long-standing and widely-held belief that there is a strong 
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connection between black hole spin and relativistic jets (IBlandford fc Znajeklll977l . hereafter 
BZ). If relativistic jets are powered by black hole spin, then theory pre dicts that jet power 
will increase dramatically with increasing a^^ (ITchekhovskoy et al.ll2010l ). For low spins, the 
black hole contributes very little power; i n fact, for g^, < .4, the accretion disk apparently 
provides more power than the black hole ( iMcKinneyl l2005l ) . 



Given the low spins of XTE J1550-564 and A0620-00, it would appear that their episodic 
jets are driven largely by their accretion disks. One well-known cand idate mechanism is 
the ce ntrifugally driven outflow of matter from a disk described by IBlandford fc Payne 
( 119821 . h ereafter BP). A usefu l comparison of the operational regimes of BP and BZ is 
given by iGarofalo et al.l (|2010[ ). They show that BP is always viable, but that BZ is a more 
likely mechanism for the most rapidly rotating sources, such as the extreme-K err black hole 
GRS 1915+105 (see Sec. 6.2). In closing, we note that a statistical study by iFender et al. 



(I2OIOI ). which is based on data of uneven quality, found no evidence that black hole spin 
powers jets. 



The Fe Ka Reflection Method 



In the Fe Ka method, one determines rigco by modeling the profile of the broad and 
skewed iron line, which is formed in the inne r disk by Doppler effects, light bending, rel- 



ativistic be aming, and gravitational redshift (IFabian et al.l I2OOOI : iReynolds fc NowakI 12003 



Millerll2007l ). Of central importance is the effect of the gravitational redshift on the red wing 
of the line. This wing extends to very low energies for a rapidly rotating black hole {a^, ~ 1) 
because in this case gas can orbit near the event horizon, deep in the potential well of the 
black hole. Relative to the CF method, measuring the extent of the red wing in order to 
infer a^, is hindered by the relative faintness of the signal. However, the Fe Ka method has 
the virtues that it is independent of M and D, while the blue wing of the line even allows an 
estimate if i. As noted in Section 1, this method, while applicable to both classes of black 
holes, is presently the only viable approach to measuring the spins of supermassive black 
holes. 

For stellar-mass black holes, in addition to the thermal disk component of emission 
(which is central to the CF method), a higher-energy power-law component of emission is al- 
ways observed (e.g., see Fig. 6). For supermassive black holes in active galactic nuclei (AGN), 
this power- law component is dominant, and it is thought to be produce d by inverse Comp- 
ton scattering of soft th ermal photons in a hot [kT ~ 100 keV) corona (IReynolds fc Nowak 
2OO3I : iDone et al.ll2007| ). Meanwhile, the disks of both stellar-mass and supermassive black 
holes are too cool [kT ~ 1 and 0.01 keV, respectively) to produce the observed Fe Ka 



- 14 - 



emission line. Rather, this dominant hne, and a host of other hnes, are generated via X-ray 
fluorescence as a result of irradiation of the disk by the hard, coronal power-law component. 
The complex spectrum so generated is referred to as a "reflection spectrum." In order to 
determ ine the spin using the Fe- line method, one must model the reflection spectrum in 
detail JRoss fc Fabian! boosl . [2OO7I ). 

The spins of several stellar-mass black holes have been measured using the Fe Ka 



metho d. An early suggestion of high spins for two black holes was made by iMiller et al. 



J2OO2I. l2004j) and p reliminary results for a total of eight stellar-mass black holes are given in 



Miller et al. 



(120091). Ot her important pape rs on the spin s of stellar-mass black holes include 



Reis et al.l ( l2008l . |2009[ ). For a review, see iMillerl ( 120071) . Very recently, we teamed up with 



Fe Ka experts to measure the spin of XTE J1550-564 (jSteiner et al.ll2010b( ). The spin esti- 
mate obtained using the Fe Ka method is a* 
CF value (see Table 1). 



0-5^-a22' which is quite consistent with the 



The spins of se veral supermassive black holes have been reported, which range from 



^ 0.6 to > 0.98 (IBrenneman fc Revnoldd I2OO6I : iFabian et al.l 120091 : ISchmoll et al.l 12009 



Miniutti et al.ll2009l ). By fa r, the most well studied o f these is the Seyfert 1 galaxy MCG-6 



30-15 (for background, see iRevnolds fc Nowak 



2003). The 6.4 keV Fe Ka hne of this AGN 



is extremely broad and skewed. IBrenneman &: Reynolds! (|2006l ) and iMiniutti et al.l (|2009[ ) 
show that the red wing extends downward to below 4 keV and conclude that a* > 0.98. 



8. Conclusion 

We have discussed the only two established classes of black holes, stellar-mass and super- 
massiv^, and the two principal approaches to measuring their spins, the continuum-fitting 
and Fe-Ka methods. Spin measurements for eight stellar-mass black holes are presented, 
and these data are used to argue that the high spin of M33 X-7 is natal, and that at least 
some relativistic black-hole jets are powered by their accretion disks, not the spin of the 
black hole. 

Two aspects of this work excite us greatly. 

First, by measuring a black hole's spin, after earlier measuring its mass, we are able to 
completely characterize the intrinsic properties of each of the black holes we study. The No 



^There is evidence for a class of intermediate-mass (100 — 10^) black holes ([Miller fc Colbert! 120041 ) . 
However, to date there existence remains uncertain because no direct and confirming measurement of mass 
has been made. 
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Hair Theorem states that a macroscopic black hole, regardless of how massive it may be, 
is described by just two parameters: M and ajf|. But is the No Hair Theorem really true? 
The only way we will answer this question is by first measuring M and a^, for a good sample 
of black holes, and then testing whether the Kerr metric corresponding to these values of 
M and a^, is completely consistent with all observables that are sensitive to the space-time 
near the black hole. In a sense, we have attempted the first test of the No Hair Theorem by 
measuring the spin of XTE J1550-564 by two independent methods, the continuum fitting 
method and the Fe Ka method (Sec. 7), and finding agreement. However, the errors in the 
two measurements are currently rather large, and we do not yet understand all systematic 
sources of error, so it would be premature to claim success. But this example provides a 
taste of how astrophysics can contribute to deep questions in physics. 

The other aspect that excites us is all the areas of astrophysics that our work ties to, 
e.g., the connections that are beginning to be made between measurements of spin and the 
phenomenology and theory of relativistic jets (Sec. 6.3), and the processes that lead to black 
hole formation (Sec. 6.1). We hope to see spin data used to help constrain models of gamma- 
ray bursts, black hole formation, black-hole binary evolution, high- and low-frequency X-ray 
oscillations, black hole X-ray states and state transitions, models of X-ray coronae, etc. 

These two strong motivations stimulate us to continue firming up the measurements of 
black hole spin, with the goal of amassing a good sample of a total of 12-15 measurements 
during the next several years. 

We conclude by noting that it is reasonable to expect LIGO, LISA and other gravitational- 
wave observatories to provide us with intimate knowledge concerning black holes. However, 
these breakthroughs are still some years in the future whereas astrophysical techniques are 
providing information on black holes today. Also, gravitational-wave facilities are unlikely to 
help us understand MHD accretion flows in strong fields, or the origin of relativistic jets, or 
the formation of relativistically-broadened Fe lines and high-frequency quasi-periodic oscilla- 
tions, etc., phenomena that are now routinely observed for black holes. In short, observations 
of accreting black holes show us uniquely how a black hole interacts with its environment. 

There is no straight path to unlocking the mysteries of black holes, probing the extreme 
physical conditions they generate, and understanding their importance to astrophysics and 
cosmology. It behooves us to explore widely because it is the synergistic exploration of 
all paths that will enlighten us. Therefore, it is important to maintain balance between 
gravitational-wave and electromagnetic studies of black holes. 



^In principle it could also have an electric charge, but astrophysical black holes are unlikely to have 
enough charge to be dynamically important. 
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Table 1. Spin Results to Date for Eight Black Holes' 



Source Spin a* Reference 



1 


GRS 1915+105 


> 0.98 


McClintock et al. 2006 


2 


LMC X-1 




92+°-°^ 


Gou et al. 2009 


4 


M33 X-7 




0.84 ±0.05 


Liu et al. 2008. 2010 


3 


4U 1543-47 




0.80 ±0.05 


Shafee et al. 2006 


5 


GRO J1655 


-40 


0.70 ±0.05 


Shafee et al. 2006 


6 


XTE J1550- 


-564 


U-'J'±-0.28 


Steiner et al. 2010b 


7 


LMC X-3 




< 0.3'^ 


Davis et al. 2006 


8 


A0620-00 




0.12 ±0.18 


Gou et al. 2010 



'^Errors are quoted at the 68% level of confidence. 

^'Provisional result pending improved measurements of M 
and i. 
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Fig. 1. — Scale drawings of 21 black hole 
Mercury distance (0.4 AU) are indicated at 
giant GRS 1915+105 with an orbital period 
orbital period of 0.2 days. The shapes of the 
and the black hole is located at the center 
inclination of the binary to our line of sight 
an inclination angle of z = 0° corresponds 
plane of the sky and is viewed face on (e.g., 
J1819.3-2525). 



binaries. The size of the Sun and the Sun- 
the top. The systems range in size from the 
of 30.8 days to tiny XTE Jl 118+480 with an 
tidally distorted stars are accurately rendered, 
of the accretion disk (see key in inset). The 
is indicated by the tilt of the accretion disk; 
to a system whose accretion disk lies in the 
i = 21° for 4U 1543-47 and i = 75° for SAX 
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Fig. 2. — Radius of the ISCO in units of GMj(? versus the black hole spin parameter. 
Negative values of correspond to retrograde motion, with the black hole spinning in the 
opposite sense of the disk. Stellar black holes are expected to have prograde spins (a* > 0) as 
a consequence of their formation in a binary system, whereas the spins of supermassive black 
holes, which are condition ed by galaxy merger events, may be either prograde or retrograde 



le.g- 



Garofalo et al.lboioh . 




Fig. 3. — [top) Accretion-disk luminosity in Eddington-scaled units ( for M = 10 Mc ^-, ) versu s 
time for all the 766 spectra considered in the study of LMC X-3 by ISteiner et al.l (j2010al ). 
The downward arrows show RXTE data which are off scale. Data in the unshaded region 
satisfy our thin-disk selection criterion L/L^^d < 0.3 (Sec. 3). The dotted line indicates the 
lower luminosity threshold (5% L/LEdd) set to avoid confusion with strongly Comptonized 
data, {bottom) Fitted values of the inner disk radius r^n = Ri^/ {GM / c^) are shown for 
thin-disk data in the top panel that meet the selection criteria of the study (a total of 411 
spectra). Despite large variations in luminosity, r^^ remains constant to within a few percent 
over time. The median value for just the 391 selected RXTE spectra is shown as a red dashed 
line. 
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Fig. 4. — Luminosity profile s from GRMHD simulations (solid lines) compared with those 
from the iNovikov fc Thornd model (dashed lines) for a* = 0, 0.7, 0.9 and 0.98 (bottom to 
top). The ISCO is located at the radius where the NT disk luminosity goes to zero. 
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Fig. 5. — (a) Allowed values of black hole mass and distance for GRS 1915+105 fall within 
the shaded triangular region (see text). (6) Six estimates of the distance to GRS 1915+105 
are shown. They range from below 7 to above 12 kpc. We are working toward a 10% 
trigonometric distance. Two hypothetical and possible outcomes of our VLBA observations, 
labeled VLBAl and VLBA 2, are indicated at the to p of panel a. For references on distance 
estimates, see Figure 18 in iMcClintock et al.l ( l2006l ). 
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Fig. 6. — Model fits for a pair of spectra of XTE J1550-564. (left) A spectrum with a 
strongly dominant thermal component; shown are the data, the fit to the data and the fitted 
thermal component, (right) A strongly Comptonized spectrum. Note the intensity of the 
power- law componen t relatiy e to its intensity in the left panel. For details, see Figure 4 and 
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